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Abstract 



We calculate the spatial correlation of electrons and positrons emitted by in- 
ternal pair conversion of Coulomb excited nuclei in heavy ion collisions. The 
alignment or polarization of the nucleus results in an anisotropic emission of 
the electron-positron pairs which is closely related to the anisotropic emission 
of 7-rays. However, the angular correlation in the case of internal pair con- 
version exhibits diverse patterns. This might be relevant when investigating 
atomic processes in heavy-ion collisions performed at the Coulomb barrier. 
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I. INTRODUCTION 



Heavy-ion collisions at energies in the vicinity of the nuclear Coulomb barrier lead to an 
alignment of the colliding nuclei. This implies that the magnetic substates are no longer 
equally populated. To describe deexcitation processes following heavy-ion collisions such 
as 7-ray emission or internal conversion, we have to account for this specific population by 
weighting the transition matrix elements with the occupation probability of, rather than 
just averaging over the decaying substates. 

The population of the various nuclear substates is incorporated in the formalism by introduc- 
ing the density matrix of the excited quantum system or, in the case of rotational symmetry 
of the problem, by a set of statistical tensors which obey the same transformation law as 
the spherical harmonics. This concept enables us to treat the polarization or alignment of 
excited nuclei appropriately. First calculations of the angular correlation of electrons and 
positrons emitted in internal pair conversion taking into account the alignment of nuclei were 
accomplished by Goldring, Rose and Warburton These calculations were performed 

within Born approximation, neglecting the influence of the nuclear charge on the outgoing 
electron and positron. But for internal pair conversion (IPC) of highly charged nuclei the 
Born approximation is not justified as can be verified by the corresponding positron spectra 

II- 

Therefore we reconsider in the following the internal pair conversion of heavy nuclei which 
are aligned, e.g., by Coulomb excitation or transfer reactions. We determine the angular 
correlation of the emitted electron and positron with respect to a reference axis in space. 
As already known for the angular correlation of 7 rays, the problem will be simplified if we 
choose a coordinate system in which the density matrix is diagonal. The statistical tensors 
depend as well on the choice of the coordinate system. If the entries of the statistical tensors 
are given in a specific coordinate system, we are able to calculate the angular correlation 
with respect to the z axis of this system. 
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The occupation probabilities of the magnetic substates caused by Coulomb excitation can 
be calculated with, e.g., the COULEX code of K. Alder and A. Winther 0. However, one 
should take into account the change of population by electromagnetic transitions from higher 
lying states. Special attention should be paid to a proper choice of the coordinate system 
when dealing with the COULEX code 03- For pure Coulomb excitation we will assume 
the z axis to point along the asymptotic target recoil axis. With respect to this axis the 
excited nuclei may exhibit prolate or oblate alignment. 



II. DENSITY MATRIX AND STATISTICAL TENSORS 



The density matrix - and for spherical symmetry the set of statistical tensors - is the appro- 
priate tool for including statistical properties such as occupation probabilities of quantum 
mechanical states into the calculations. Here we briefly summarize the essential properties 
of the density matrix and subsequently turn to the concept of the statistical tensors, which 
obey the same transformation law as irreducible tensors. For the density matrix pMiM'.{Ji) 
of dimension (2Jj + 1) x (2Jj + 1) we note: 

1. The density matrix is hermitian: 

Pm'.mAJi) — PMiM[{Ji) , 

1 i I 

2. The trace of the density matrix equals one: 

Tr{p} = 1 $>W4) = 1 , 

Mi 

3. Tr{p 2 } < 1 , and Tr{p 2 } = 1 4$ the system is in a pure state. 

We define the statistical tensors pfi as irreducible tensors of rank n with v = —n, . . . , n: 

pMiM'XJi) ■ (1) 



Mi,M' 



y Mi M[ -V j 
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The argument Jj reminds us that n is related to the angular momentum of the magnetic 
substates by < n < 2J$. From the normalization of the density matrix it follows po(Ji) = 



i/v/2J7+T. 

The density matrix has (2Jj + l) 2 independent components. To describe a system by sta- 
tistical tensors instead of the density matrix, we need 2 Jj + 1 density tensors of rank n = 
up to rank n = 2 Jj. Since the density tensor of rank n has 2n + 1 components we get again 
En=o (E"=_„ l) = (2 J + l) 2 independent components. 

The statistical tensors transform under rotations according to 

p [ : ] (Ji) = Y.v [ X*(*)p [ M) (2) 

where the Wigner rotation matrix of rank n is denoted by T>^ and the set of Euler angles 
by a. In defining the Euler angles we follow Rose || and Eisenberg&Greiner [|IIJ]. For 



systems with rotational symmetry it is thus more advantageous to employ the concept of 
the statistical tensors when incorporating statistical statements concerning the system. The 
components of the statistical tensors are changed under rotations and so are the occupation 
numbers of the magnetic substates. 

From the set of 2Jj + 1 statistical tensors we obtain the density matrix by utilizing the 
relation: 



PM iM >{Ji) = (-i) J ^ M >E^TT 

n.u 



y M[ -Mi -V j 



P [ u\^) ■ (3) 



For certain symmetries of the system we can reduce the independent components of the 
statistical tensors. Here we list the consequences for the statistical tensors in three special 
cases which will become relevant for us: 



1. In the case of axial symmetry the density matrix is diagonal, its diagonal components 
are just the probabilities for the occupation of the corresponding magnetic substates 
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PM t Mi — PMi- 



-Mi j 

One can always choose a basis such that the density matrix is diagonal, but in general 
this will not be a basis of wave functions with good angular momentum. 

2. In the case of spherical symmetry, there is no direction singled out in space. The 
density matrix is proportional to the identity matrix. The diagonal elements are given 
by PMiMi = 1/(2 Jj + 1). All statistical tensors vanish with exception of the tensor of 
rank 0, i.e., 



yJUi + 1 

3. From Eq. (|I|) it can be shown that for alignment of the nuclear states, defined by 
PMi — P-Mi, the statistical tensors of odd rank vanish. 

III. ANGULAR CORRELATION OF 7-RAYS 

Before we enter into the calculations concerning the angular correlation in internal pair con- 
version we summarize some results already known for in-beam 7-ray spectroscopy. This will 
help us to interpret the angular correlation pattern in the case of internal pair conversion. 
The angular correlation of photons emitted after Coulomb excitation is essentially deter- 
mined by the statistical tensors, i.e., by the occupation numbers of the magnetic substates 
of the decaying nucleus. In choosing a reference axis for which the density matrix is diagonal, 
just the 0th components of all statistical tensors survive and we obtain for the transition 
probability the well-known relations: 
dP 7 2acu 



dVL ^2Ji + 1 



Vj T \ L )\* Yl Fi{ LLJ f J i)Po ] {Ji) p i(cos'd) , (4) 

1 even 



for a transition of parity r = E/M and multipolarity L. V^ T \L) denotes the correspond- 
ing reduced matrix element for the nuclear transition. Here we employed the correlation 



coefficients [[T]]-|13 



Fi(LL Jf Ji) = (-) J / +J *- 1 v / 27TT V / 2J i + l 
x (2L + 1 



U L l\ 




L L I 






! 




(5) 


v 1 - 1 °) 




Ji Ji Jf 





This results in an anisotropic emission of photons with respect to the alignment axis. The 
number of minima of the angular distribution corresponds to the multipolarity of the nuclear 
transition. 



In the case of spherical symmetry, the photon emission is isotropic 



da 



2auj 



dtt 2Ji + 1 
or integrated over the solid angle Q: 



(6) 



2X + 1 



Vi T \L) 



(7) 



IV. TRANSITION PROBABILITIES FOR INTERNAL PAIR 

CONVERSION 

We turn now to the formulation of the triple correlation of the electron and positron direction 
with reference to a symmetry axis, which is taken as quantization axis. 

For a statistical ensemble of nuclei we write the transition probability for internal pair 
conversion, 

p e+e _ = 2tt I d3 p [ d V S(oj -W'-W) 

Mi,M(,M f ,X,X' 

XlJpYPM.M'U*! . (8) 

where the density matrix Pm 1 m[ represents the occupation of the magnetic substates \JiMi). 
Here we assumed a nuclear transition from a initial state | JjMj) to the final state \JfMf) 
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where the initial state is populated according to the density matrix puihii- Since we do not 
require the density matrix to be diagonal the summation extends over both, Mj and M-. 
The 5 function ensures energy conservation: The transition energy u is transfered to the 
electron (total energy W) and to the positron (total energy W). The summation is taken 
over the spins and the momenta of the outgoing leptons. 

The matrix element for internal pair conversion is written in lowest order of a in the retarded 
form: 

U P \ = -at J dV n J dV c (p n (r n ) p e (r e ) - j n (f n ) ■ j e (r c )j j - , (9) 



r c being the electronic coordinate, r n the nuclear coordinate. 

Since we neglect in our work the penetration of the electron wave functions we do not have 
to specify the nuclear transition charge and current densities p n and j n . The electronic 
transition charge and current densities read 

p e = Ipi , Je = 4>j a Ipi (10) 

(a is the 3- vector of the spatial Dirac matrices in the standard representation). These 
expressions are evaluated utilizing the scattering solutions, see Eqs. ( |/^5| , |/^8|) in the appendix, 
for the electron and positron wave function in order to define the emission direction and thus 
an opening angle. Inserting the spherical wave expansion of these wave functions results in 
a decomposition of the matrix element, Eq. (Q), 

^=EE°Sy ^ ^ • (11) 

Uk'h'kh denotes the transition matrix element which has the same structure as U p \, but is 
evaluated using the spherical spinor solutions of the Dirac equation, Eqs. (Al). This matrix 
element was calculated in IB. Here we cite the result: 



U K ,„, K „ = iTriau{-l) J f- M f 



—Mf M Mi 



x(-iy 



i 



-// M /j, 



V^ T \L) is just the reduced nuclear matrix element of Eq. (|7p and 



AnJL(L + 1) 



x 



-1/2 1/2 



with the parity selection rule 



/ + /' + L + A(r) = mod 2 , 



A = for r = el 
A = 1 for r = maern 



(12) 



(13) 



(14) 



R$ K contains the integration over the radial electron wave functions and will be defined 
later. 

Inserting this matrix element into the pair conversion probability, Eq. (|8]), yields 
p e+e _ = 2tt J2 Pm]m> J2 [ dWdn I dW ' dn 6 ( u -W'-W) 



Mi,MiM f A, A 

x X] X! -^k/w Uk'h'kh Uftp, 



where we abbreviated 



(15) 



(16) 



and 



(17) 



From Eq. ([TJJ) we obtain the differential pair conversion probability with respect to the 
kinetic positron energy E = W — m and the solid angles of both electron, Q', and positron, 

a 
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u-2m 



dE dQ dQ 



dE dQ dQ' 



(18) 



The integration over the electron energy W is trivially performed because of the 6 function. 
From this relation it is obvious that we may proceed from the solid angles Q to Q by choosing 
another reference axis in space. The integrand in Eq. ( jl8|) is invariant under rotations since 
the Jacobian of this transformation equals 1. The integrand should thus be represented by 
a series of triple correlation functions which are defined in Eq. (|22]) . 



Inserting the explicit expressions of the coefficients A and B, Eqs. 
following expression for the differential pair conversion probability 

d 3 P fi - 



dEdQ dQ 
( 

x 



7 = 8(^) 2 |^)(L)| 2 E £ EpSm; 

M,M Mi,M'- M f 



Jf L Ji 



\ 



( 



y-M f M Mi 



Jf L J4 



i —Mf M Ml 



x 



E (-!) /+? M ^( L ) Mit(L)pf + 1 VV + 1 



xy^j + 1 V2J + 1 

x exp(i[(5'(W, /«') - 8'(W, R') + 5(-W, k) - 5(-W, re)]) 



x 53 v / 27 7 TIv / 27ny/w(ft P ')}/«(^) 

r,I,a,a' 

( - - r \ ( 



X 



1/2 -1/2 



1/2 -1/2 



E (-ir* ,+l 



; i ^ 



v 



X 



V 



f L j 
-fjf M \i 



Ji' M Ji 
M ?' i' I' \\ 3 3 I 



3 3 



leads to the 



Here we inserted Eq. fll2"|). Introducing the statistical tensors we are left with 

d 3 P P - 



dE dQ dQ 



7 = 2(4™) 2 |V;W(L)| 2 (-1)^ 
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Ji+L+1 



(19) 



x£V2^(-irpM(j,; 

J i J i J f 




x ^ v / 2/ + 1 y/2I' + 1(-1) / ' 

/ 

a,a' \ a a' —v 



I I' n 



x ^ s]\kk'RW\ 

;/ 1 ttiI t\ r/ / ttt/ - / 



I I' n 



}M£l(L)M£>*(L) 



(r)*. 



x exp(i[5'(£', «') - 6'(E', k) + 5(-E, k) - 6(-E, re)]) 



X 



1/2 -1/2 



1/2 -1/2 



(20) 



This is the most general form for the pair conversion probability. Now we assume that 
we are dealing with internal pair conversion of aligned nuclei iv = 0). We may choose an 
appropriate coordinate system by transformation of the spherical harmonics: 



Y Ia (a p 



£ exp(i«0) exp(i/35) Y^(0, 0) 



exp(ia0) <4fo($) 



V 4vr 

]T exp(-ia0) d [ %(>&) exp(i/3'5) Y r p,(Q, 0) 



(21) 



Here G denotes the opening angle of the electron-positron pair, $ is the polar angle of the 
positron with respect to the symmetry axis, and the dihedral angle 5 indicates the rotation 
of the electron-positron plane around the positron axis (Fig. 2b). Please note, that the 
convention of |]^§ differs in the definition of the angles from the one employed here. 



This enables us to define the triple correlation function by 



p II , n (#,e,5) = y £ 



( I V n\ 



a —a 



Yi a (Qp) Y//_ a /(f2 p ') 
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y / 27TTy / 2^ + 1 
4tt 



1 ' I V 



0/3/?' 



d% (0) 43(9) exp(i/3'5) . 



(22) 



Our triple correlation function is related to the one introduced by Biedenharn [|Tl|] by a 
factor 4vri- / - / '- r V((2J + 1)(2/' + l)) 1 / 2 . 

The pair conversion probability is normalized by the probability for 7 emission, Eq. (^), 
which yields the pair conversion coefficient: 



A 2 - = 2 r' r 2L + 1) (2j, + 1) ( -D^-.«« 



d£dcos9dcos0d5 L(L + 1) 

L L n 



X 



£(2/+l)(2/' + l)(-l/' 



xE(-if 

/3' 



/?' -/?' 



dfc ] (tf) rf [/ ' ] /?'0(e) exp(i/3'5) 



x 



2 (-1) s+3 '\k k'kk' I 



(t)*. 



J y £ 

I V n 

x exp(i[<f'(W", K ') - 5'(W, R') + 5(-W, k) - 5(-W, «)]) 



x 



1 f J 



A( 



x 



1/2 -1/2 

j f L 
1/2 -1/2 



j J I 



\ 



( , ,/ r \ / 1 * r\ 



1/2 -1/2 

3 f L 
1/2 -1/2 



(23) 



Here we inserted the explicit expressions for the electronic matrix elements, Eq. (|lq) . Inte- 
gration over the azimuthal angle is trivially performed resulting in an additional factor of 
2tt. 

The radial matrix elements read for electric pair conversion (parity (— ) L ) 



12 



R ( X = L(Ri + R2 + R3- Ri) + (« - k')(R 3 + Rt) (24) 

and for magnetic pair conversion (parity (— ) L+1 ) 

R^^^ + k'KRs + R*) . (25) 

The radial integrals introduced in these equations are taken over products of the radial 

electron wave functions (|A4|) and the Hankel functions of first kind, h^\ur): 

00 

Ri= J &rr 2 g w , tK ,{r)g_ w , K (r)h { l\ujr) , 

»..;,.,,,„„-..,,«"„. 



00 



00 



R \ 





00 

4 = /*^^)/.^)^) 



00 



Re = Jdrr 2 f wltKl (r)g^ K (r)h { l\ujr) . (26) 


In the case of a point-like nucleus these integrals can be rewritten in terms of F2 func- 
tions which can be evaluated numerically. For the representation of the nucleus as a 
homogeneously charged sphere the radial integrals are computed using a Gauss-Chebyshev 
quadrature |14|. The Whittaker functions which occur in the expressions for the electron 
wave functions are computed with the COULCC code of [I5|. Since the integrands are oscil- 



lating functions it is advantageous to deform the integration contour in the complex plane 
in such a way that it runs along the imaginary axis ||. Since the electron wave functions 
have the asymptotic behaviour exp(fpr) while the Hankel functions behave like exp(kjr), 
where to = W + W, the integrand for large r assumes the form: 

exp(i[-p-p' + W + W'}r) . (27) 
13 



For r complex with the imaginary part going to infinity, our procedure thus guarantees that 
the integrand falls off quite fast. In most cases at r = 20000 fm the integrand is smaller 
than 10 -5 of its maximum value. 

We want to consider the angular correlation for two special cases: 



I. If we integrate Eq. ( |2"3"| ) over the positron polar angle $ and the dihedral angle 5 the 
remaining function depends only on the opening angle G of the electron-positron pair. In 
this case only the n = contribution survives. We get the opening angle distribution as a 
series of Legendre polynomials which was already calculated in [ |16[] : 

d 2 {3 



dE d cos 6 



a/ Pi (cos 0) 



(28) 



The expansion coefficients are given by 



a/ = L^TrT ( - )L+m(2/ + 1) £ \^'KR'\R^ K {L)Rti 

\ / k,k',k,k' 

x exp(i[8'(W, k') - 5'(W, k') + 5(-W, k) - 6(-W, R)]) 

/..-.. T \( . - A 



X 




(29) 



They have to be evaluated numerically. The same result is achieved if one assumes that the 
initial nuclear substates are equally populated. 

At this point we apologize for giving an incorrect expression for the opening angle distribu- 



tion in [IB| which was caused by employing the wrong set of scattering solutions. This error 
resulted in the wrong sign of the scattering phase shifts of the positron. The opening angle 
distribution showed the right qualitative behaviour but wrong conversion probabilities. The 
statement, that the maximum of the distribution shifts from 0° to 180°, if one considers 
overcritical nuclear charges [Z > 173), remains unchanged. This error appeared also in 
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the expression for the electric monopole (EO) conversion. One should reverse the sign of 
the scattering phase shifts of the positron. The pair conversion coefficient for the electric 
monopole conversion reads: 

dV -^(1 + ecosG) (30) 



d£dcos6 2 dE 

where drj/dE is the differential pair conversion coefficient || — which remains unchanged- 
and e is the corrected anisotropy coefficient: 



e = 2 c\ + c\ 1 C0S ( A + 1 - 1 ) ( 31 ) 



with 

A KK , = 5'(W',k)-5'(W',k') + 5(-W,k) -5(-W,k') . (32) 



C + i, C-i are defined by 



a 



' .. f-wA 7 *) f\V>A r ) r ^ n 

hm -— ■ — ■ tor k > 

g-W,K( r )fw',n( r ) f ^ n 

hm r— ; — ■ tor k < 



(33) 



where / (/') and g (g f ) are the radial wave functions of the Dirac spinor of the positron 
(electron). In numerical calculations these constants are evaluated at the nuclear radius. 5 



and 5' are the corresponding Coulomb phase shifts for an extended nucleus [17 



II. If we integrate Eq. ( |2"3"D over the opening angle G of the electron-positron pair and the 
dihedral angle 5 we end up with 

d 2 {3 



dE d cos i? 



^&„P n (cost9) 



(34) 



where the coefficients read 

K = W PL + 1) ( _ )J ,_ J , + , V ^ TT | L L n 

I "^i J% Jf 
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L L n 



3 3 3 



x exp(i[<?(-W, k) - 6{-W, re)]) 



3 3 n 
1/2 -1/2 



3 3' L 
1/2 -1/2 



3 f L 
1/2 -1/2 



(35) 



This corresponds to the experimental setup where one is just interested in the angular 
distribution of the positron emitted in internal pair conversion of an aligned or polarized 
nucleus. 



V. RESULTS 

In the following we will discuss the characteristic properties of IPC angular distributions 
using a few representative results. Since we are interested in Coulomb effects they all refer 
to a uranium-like nucleus (Z = 92). The chosen energies and multipolarities are generic and 
are not intended to represent particular nuclear transitions known from experiment. 

The opening angle distribution of electron and positron emitted by IPC depicts for electric 
transitions the typical pattern: it has its maximum at O = 0° and its minimum for = 180°. 
For magnetic transitions in heavy nuclei, however, the situation might be different. Fig. la 
and lb depict the opening angle distribution for an El and a Ml transition of a uranium- 
like nucleus as a result of our distorted wave Born approximation (DWBA) in comparison 
with the Born approximation (BA). This demonstrates how the angular correlation of the 
electron-positron pairs is influenced by the strong Coulomb field of the nucleus. In Fig. lb 
we plotted also the opening-angle distribution for the Ml transition taking into account 
the finite extension of the nucleus under consideration. This verifies that the magnetic 
transitions — and especially the Ml transition — are very sensitive to the charge distribution 
of the nucleus For the El transition in Fig. la, on the other hand, the effect of the finite 
nuclear size amounts to less than 0.1 %. 
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In the following we discuss the triple angular correlation of electron and positron for IPC 
of aligned nuclei. We take the symmetry axis as quantization axis as in Eq. (|23|). The 
opening angle B of electron and positron, the polar angle $ of the positron and the dihedral 
angle 5 form a complete set of angles to fix the emission directions of electron and positron 
with respect to the symmetry axis. The angles which describe the directions of the emitted 
leptons are displayed in Fig. 2b. Note that our choice of the coordinate system is different 
from that introduced in the Born approximation calculations of fl|-|3[ in which $ denotes 
the polar angle of the intermediate photon. However, since the Coulomb field disturbs the 
momentum balance we cannot determine the momentum of the intermediate photon from 
the momenta of the outgoing leptons, which would be necessary to calculate the photon 
polar angle. 

Depending on the experimental setup and reactions various coordinate systems may be 
established in which the statistical tensors are determined. Here we concentrate on the 
Coulomb excitation of heavy ions in collisions with beam energies at or below the Coulomb 
barrier. In this case one usually chooses a coordinate system, where the z axis is pointing 
along the apex line of the scattering hyperbola towards the projectile and the x axis is 
perpendicular to the scattering plane (Fig. 2a). The y axis is then chosen such that the y 
component of the projectile velocity is positive In the sudden approximation it can 

be shown that the nuclear states are excited with a population of the magnetic substates 
reaching a maximum around Mj = 0, i.e., the nucleus is aligned in the plane perpendicular 
to the z axis (asymptotic recoil direction of the target) |13|]. This is called oblate alignment. 
Taking into account the de-excitation of the nucleus by 7 cascades starting from high-spin 
the oblate alignment changes into a prolate alignment with respect to the z axis for the 
low-spin states. 

If the collision energy is increased the nuclear alignment changes to a polarization with 
respect to a reference axis perpendicular to the scattering plane Classically this 

corresponds to the situation where the drag caused by surface friction puts the nuclei into 
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a spinning motion. 

After having chosen a coordinate system and having determined the degree of alignment or 
polarization for the Coulomb excited nuclei — the corresponding statistical tensors can be 
calculated with, e.g., the COULEX code of — one can employ Eq. (|23D to determine the 
angular distribution of the electron-positron pairs emitted by internal pair conversion of these 
nuclei. We plot in Fig. 3 the spatial correlation of the electron-positron pairs with respect to 
the reference axis assuming oblate alignment of a uranium-like nucleus. From the spectrum 
of the emitted pairs |5,19] we know that for large- Z nuclei the pair emission probability 



increases towards the maximum positron energy. Thus the angular correlations are plotted 
for a case where nearly the full transition energy (minus the electron rest mass) is transferred 
to the positron. One recognizes a strong dependence of the pair conversion probability on 
the polar angle of the positron with respect to the reference axis. This behaviour resembles 
the anisotropic emission of the intermediary photon [f[3f| . The angular distribution depends 
weakly on the dihedral angle 5 of the electron-positron pair (Fig. 4). For transitions between 
nuclear states of high angular momentum the opening angle distribution does not change 
drastically when the positron polar angle is varied. 

In order to elucidate the influence of the statistical tensors, i.e., the occupation of the initial 
nuclear state on the angular correlation of the emitted electron-positron pairs, we present 
the angular distribution with respect to the polar angle of the emitted positron. Fig. 5 
shows the polar angle distribution assuming El, E2 and E3 transitions to the + ground 
state of nuclei which exhibit oblate alignment. Fig. 6 displays the polar angle distribution 
for a El transition to the + ground state of a nucleus for oblate and prolate alignment and 
for polarization. 
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VI. CONCLUSION 



Angular correlations are very sensitive to the underlying process. They may reveal a plethora 
of signatures for nuclear transition which allow for an identification of the transition as well 
as for the study of the properties of excited nuclei. Especially for large- Z nuclei one cannot 
rely on the validity of the Born approximation which becomes exact in the limit Z — > 0. 
One rather has to perform the calculations with the relativistic scattering wave functions 
for both electron and positron. These wave functions take the Coulomb distortion caused 
by the nuclear charge into account. 

For magnetic transitions of large- Z nuclei one has also to account for the finite-size effects. 
Especially Ml transitions are very sensitive on the extension of the nucleus. In order to 
study magnetic IPC we approximated the nucleus undergoing the transition by a homoge- 
neously charged sphere. The angular correlation of electron-positron pairs with respect to 
a given axis in space depend on the statistical tensors which reflect the population of the 
nuclear magnetic substates. The conversion probability changes drastically when either the 
opening angle of the pair or the polar angle of the positron is varied while the dependence on 
the dihedral angle is rather weak. Our calculations allow to make quantitative predictions of 
this behaviour which qualitatively might have been anticipated from the 7-ray spectroscopy 
performed in heavy-ion collisions. Furthermore, the measurement of the spatial correlation 
of electron-positron pairs can be employed to obtain additional information about the nu- 
clear transition. E.g., not only the multipolarity but also the parity of the nuclear transition 
can be measured in this way. 



Acknowledgement: This work has been supported by the BMBF, by the Deutsche For- 
schungsgemeinschaft (DFG), by GSI (Darmstadt), and by the REHE programme of the 
European Science Foundation (ESF). 



19 



APPENDIX A: ELECTRON WAVE FUNCTIONS FOR POINT-LIKE AND 

EXTENDED NUCLEI 



In our calculations we employed the following form of the spherical continuum wave functions 



of the electron moving in the Coulomb field of a point-like nucleus |20 

9w, K ( r ) Xw(r) 1 

Xw,KA r ) 



gw, K [ r )Xw[r 



The spinor spherical harmonics are defined as 



m A 



J 



where the basis spinors are given as usual by 
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Defining the relativistic Sommerfeld parameter as y — —ZaW/p where p = \/W 2 — m 2 , the 
radial wave functions read for a point nucleus 



9wA r ) 



fwA r ) 
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with 
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For an extended nucleus we construct the continuum solutions as in |I7| by employing a 
power series ansatz for the electron wave function inside the nucleus which is matched to 
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the linear combination of wave functions to the Coulomb potential at the nuclear radius. 
From the matching condition the normalization factor and the phase shift can be deduced. 

The continuum solutions of the Dirac equation can be written as wave functions which 
asymptotically represent plane waves of momentum p and spin A. These wave functions are 



obtained as a series expansion into spherical harmonics pOjpT 
For positive energies this expansion reads: 



(±) 

W,p,X L^i^KiX A.W,K,H 
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with the coefficients 
,(±) - 



1 



a) 



and the Coulomb phase shift 



and for negative energies (— W < 0) |fH| 



with the coefficients 
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fe(±) = 



/(-«) \ 



\ 



m A 



and the Coulomb phase shift 
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S(-W,k) = r] - argT(7 - iy) - -7 . 



These wave functions obey the normalization condition 
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In the case of an extended nucleus the phase shifts have to be determined numerically. 

21 



APPENDIX B: EVALUATION OF THE COEFFICIENTS A AND B 



We start with the coefficients A which contain the electron scattering wave function and its 
complex conjugate. Inserting the explicit form of these coefficients, Eq. (|A6|), we get 



= £ ex P (i[5'(W',K')-5'(W',R')]) 
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which can be transformed into 
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Additionally we get the parity selection rule 



-// /i' —a! 



I'(k') + I'(r') +I' = mod 2 



(Bl) 



(B2) 



and the angular momentum selection rule: 

\1'(k')-1'(R')\ <I'<1'(k') + T(k') . 

Next we proceed to evaluate the coefficients B, which are composed from the positron 
scattering wave functions, Eq. (|A9|) : 
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This can be rewritten as 



B K , m = (-iy +l ' 2 ^= exp(i[5(-W,K) - 5(-W,R)\) 

V47T 



x J2j + 1 V2j + 1 



xJ2V2TTlY Ia {Q 
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Furthermore we obtain the parity and angular momentum selection rules 



1(-k) + l(-R) + / = mod 2 , 



\1(-k) - l(-R)\ < I < 1(-k) + l(-R) 
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APPENDIX: FIGURE CAPTIONS: 



Figure 1: Opening-angle distribution of electron-positron pairs emitted by IPC of randomly 
oriented uranium-like nuclei a) for an El transition, b) for a Ml transition. The transition 
energy amounts in both cases to 2 MeV, the kinetic positron energy was taken to be 800 
keV. The full lines correspond to the DWBA calculations, dotted lines display the outcome 
of the Born approximation. The effect of the finite nuclear extension in the case of Ml 
transitions can be deduced from the dashed line. 

Figure 2: a) The coordinate system which is chosen to describe the alignment of the 
Coulomb excited nuclei in heavy-ion collisions, d denotes here the scattering angle of the 
projectile in the lab system, b) Definition of the angles which are used to describe the 
directions of electron and positron in space. is the opening angle of the electron-positron 
pair, i? denotes the polar angle of the positron with respect to the quantization axis, and 
the dihedral angle 5 indicates the angle about which the electron-positron emission plane is 
rotated around the positron axis. 

Figure 3: Angular correlation of electron-positron pairs in El-IPC for a transition from 
a T to a + state, assuming oblate alignment of uranium-like nuclei in the initial state. 
The conversion probability is plotted versus the opening angle of the emitted lepton pair 
for various polar angles of the positron with respect to the z axis. The transition energy 
amounts to uo = 1800 keV, the kinetic positron energy to E — 700 keV. The dihedral angle 
5 was fixed to 0°. 

Figure 4: Angular correlation of electron-positron pairs in El-IPC for the same transition 
and energies as in Fig. 3. The pair conversion probability is plotted for fixed polar angle 
i? = 90° and several values of the dihedral angle 5. 

Figure 5: Polar angle distribution assuming El, E2, and E3 transitions to the + ground 
state of a nucleus showing oblate alignment. The transition energy amounts to 1800 keV, 
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the kinetic positron energy was fixed to 700 keV. 

Figure 6: Polar angle distribution assuming an El transition to the + ground state of a 
nucleus with oblate alignment, prolate alignment, and polarization. The transition energy 
amounts to 1800 keV, the kinetic positron energy to 700 keV. 
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Figure 6 
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